arXiv:1509.07614vl [quant-ph] 25 Sep 2015 


Least-bias state estimation with incomplete unbiased measurements 


Jaroslav Rehacek,' Zdenek Hradil,' Yong Siah Teo,' Luis L. Sanchez-Soto,^’^ 

Hui Khoon Jing Hao Chai,^ and Berthold-Georg Englert^’®’’ 

' Department of Optics, Palacky University, 17. listopadu 12, 77146 Olomouc, Czech Republic 
^Departamento de Optica, Facultad de Fisica, Universidad Complutense, 28040 Madrid, Spain 
^Max-Planck-lnstitut fiir die Physik des Lichts, GUnther-Scharowsky-Strafle 1, Ban 24, 91058 Erlangen, Germany 
*Yale-NUS College, 16 College Avenue West, Singapore 138527, Singapore 
^ Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore 
^MajuLab, CNRS-UNS-NUS-NTU International Joint Unit, UMl 3654, Singapore 
''Department of Physics, National University of Singapore, 2 Science Drive 3, Singapore 117542, Singapore 

Measuring incomplete sets of mutually unbiased bases constitutes a sensible approach to the tomography of 
high-dimensional quantum systems. The unbiased nature of these bases optimizes the uncertainty hypervolume. 
However, imposing unbiasedness on the probabilities for the unmeasured bases does not generally yield the 
estimator with the largest von Neumann entropy, a popular figure of merit in this context. Furthermore, this 
imposition typically leads to mock density matrices that are not even positive definite. This provides a strong 
argument against perfunctory applications of linear estimation strategies. We propose to use instead the physical 
state estimators that maximize the Shannon entropy of the unmeasured outcomes, which quantifies our lack of 
knowledge fittingly and gives physically meaningful statistical predictions. 
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I. INTRODUCTION 

Modern quantum technologies harness characteristic fea¬ 
tures of quantum systems to gain performance that is other¬ 
wise unattainable through classical means. This progress ulti¬ 
mately relies on the ability to create, manipulate, and measure 
quantum states. All of these tasks require a step-by-step veri¬ 
fication in the experimental procedures; this is essentially the 
scope of quantum tomography m . 

Typically, a tomographic protocol attempts to infer the un¬ 
known quantum state from the distinct outcomes of a collec¬ 
tion of measurements performed on a hnite set of identical 
copies of the system. With these limited resources, the choice 
of optimal measurements and the design of efficient recon¬ 
struction algorithms are crucial. 

When the measurement outcomes form an informationally 
complete set EHSl, the data obtained contain maximal infor¬ 
mation and a unique state estimator can be inferred. Unfortu¬ 
nately, as we probe more intricate quantum systems Cl El, 
such an informationally complete set of measurements be¬ 
comes extremely difficult to implement. In addition, the com¬ 
plete knowledge of the quantum state of a system with many 
degrees of freedom is usually not needed, as one could very 
well be interested in a few parameters only, such as the fidelity 
with respect a target state or a measure of entanglement. 

We are thus inevitably led to the consideration of alterna¬ 
tive techniques. A promising class of new protocols are ex¬ 
plicitly optimized for particular kinds of states. This includes 
states with low rank (Mil, such as matrix product states 
main], or multi-scale entanglement renormalization ansatz 
states m. The specihc but pertinent example of permutation 
invariance was also examined EMU- 

In the same spirit, several algorithms for estimating quan¬ 
tum states from incomplete measurement data have been re¬ 
ported recently llT9l - Ea . Here, we revisit the problem in the 
context of mutually unbiased (MU) measurements, which are 


known to be optimal for state reconstruction ll^ if a com¬ 
plete set of bases is to be measured. At the fundamental level, 
mutually unbiased bases (MUB) are part of the mathemati¬ 
cal framework for an explicit formalism of the complemen¬ 
tarity principle. The idea emerged in the pioneering work of 
Schwinger ll24ll . and has since been integrated into the founda¬ 
tion of quantum theory: apart from their role in tomography, 
MUB are instrumental in addressing a number of enthralling 
questions ll25l . 

A first, self-suggesting, if naive, approach could be to as¬ 
sume a uniform distribution for the outcomes of the unmea¬ 
sured bases from a set of MUB, as this seems to be consistent 
with the very physical nature of MUB, which minimize the 
statistical uncertainty hypervolumes. However, as we show 
here, this is often incompatible with the positivity required by 
quantum mechanics, and even when the resulting estimator is 
physical, it is usually not the estimator with the largest von 
Neumann entropy im. 

The bases of eigenstates of complementary observables are 
called unbiased because we cannot predict at all the outcome 
of a projective measurement in one basis if the system is pre¬ 
pared in a state from another basis—all outcomes are equally 
probable, there is no preference in our prediction. This no¬ 
tion of being unbiased must not be confused with other uses 
of the adjective. In fact, there are many different meanings 
and connotations associated with “bias”: statistical bias, bias 
in a sample, cultural bias, and media bias are but a few uses 
of the word. 

In the context of state estimation, statistical bias is of some 
importance; an estimator has a statistical bias if its average 
over all thinkable compositions of a finite sample of measure¬ 
ment results deviates from the true value. This will be of no 
concern here because we shall take for granted that the mea¬ 
sured sample is so large that statistical noise in the data can be 
safely ignored. Also, we do not have to worry about a biased 
sample, a common problem when polls are taken. 
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The simple “estimators of unbiased linear inversion,” which 
we shall introduce in Sec. |IIC[ are unbiased in the sense that 
they assume equal probability for all outcomes of the unmea¬ 
sured bases from a set of MUB; they are, however, estima¬ 
tors with a statistical bias (for data from a finite sample), a 
property that cannot be avoided IZTll if one insists, as we 
do, that all permissible estimators are physical—they must be 
bona fide density operators. As noted in the preceding para¬ 
graph, the statistical bias is irrelevant in the current context, 
and it is worth recalling that, despite the negative connotation 
of “bias,” a statistical bias is not only harmless, but can be 
rather beneficial (see Jaynes’s discussion in sections 17.2 and 
17.3 of Ref 121). 

Regarding cultural bias and media bias, we are confident 
that we do not suffer from the former and hope that we will 
not be a victim of the latter. May Bias, Melampus’s brother 
and Pero’s and Iphianassa’s husband (consecutively), guard 
this work against biased judgment. 

The plan of this paper is as follows. After introducing back¬ 
ground material and notational conventions in Sec. we il¬ 
lustrate various aspects of the said self-suggesting approach 
in Sec. with particular emphasis on unphysical properties 
of the resulting estimators. We conclude that the “estimators 
of unbiased linear inversion” are often unphysical—they are 
not estimators. Then, in Sec. HYl we introduce a natural al¬ 
ternative that keeps these estimators whenever they are phys¬ 
ical and, when they are not, replaces them with physical es¬ 
timators. Rather than being unbiased about the probabilities 
of unmeasured bases, the physical estimators are minimally 
biased. In this alternative approach, we maximize the Shan¬ 
non entropy with due attention to the physical constraints and 
so minimize the bias. This entails a simple optimization al¬ 
gorithm over the state space, in which a single equation is 
iterated. In Sec. we consider alternatives to maximizing 
the Shannon entropy. Some of them minimize the bias with 
respect to another criterion, others are characterized in a dif¬ 
ferent way. Finally, we offer a summary and conclusions in 
Sec.lVH 


II. BASIC MATTERS 


of MUB; here and below, Greek indices {a, j3, ...) label the 
d+1 distinct MUB, whereas Latin indices {k, I,...) label the 
d outcomes in each basis. 

We define MU projectors as Hak = \4’ak){i^ak\- Any two 
MU projectors satisfy the trace relation 

Tr(nQ,fc 11/3^) = SafiSkl + -(1 — Safi) (2.1) 

that states their orthonormality for a = jd and their mutual un¬ 
biasedness for a ^ fi. Besides, these projectors constitute a 
complete set of projective measurements for each a, 

d-l 

J2^ak=t. (2.2) 

k=0 

To facilitate the discussion of the physics behind incom¬ 
plete MUB tomography, the number of copies of the system 
probed by the measurement apparatus is taken sufficiently 
large, so statistical fluctuations in the measurement data are 
negligible. 

The problem we are studying here is, therefore, not one 
of state estimation sensu stricto, where the estimation of the 
probabilities from observed relative frequencies is the central 
theme. Rather, we are dealing with the problem of converting 
the probabilities into a statistical operator, which requires a 
deliberate choice when the tomography is incomplete. 

B. Complete tomography with MUB 

The d'^ + d measured probabilities 

Pak = Tr(n„fc g ), (2.3) 

which obey the d + 1 constraints J^kPak = 1, fully charac¬ 
terize the density operator g of the system. First, if only one 
basis is measured—the ath basis, say—our state estimator is 

d-l 

got — ^ '^, Potk^otk • (2.4) 

fc=0 


A. Mutually unbiased bases 

The state of a d-dimensional quantum system is specified 
by a positive semidefinite, unit-trace density operator g. Care¬ 
fully note that a different symbol p is reserved for its d x d 
matrix representation, which requires d^ — 1 independent real 
parameters for its complete characterization. If a von Neu¬ 
mann maximal test is chosen to fix d — 1 of these parameters, 
then a total of d -f 1 tests is necessary to reconstruct the state. 
This strategy is optimal when the bases in which the measure¬ 
ments are carried out are “as different as possible;” that is, 
when these bases are MU 1251 . 

Throughout this paper, we take the dimension d to be a 
prime or a prime-power integer. Then, a maximal number of 
d -f 1 MUB l29l exist and can be explicitly constructed 1301- 
[39l . We denote by li/jak) the /cth ket in the ath basis of the set 


It is the most natural estimator, inasmuch as we estimate the 
probabilities of the unmeasured bases in the unbiased manner 
of pfii = Tr(n/ 3 ; Pa) = 1/d for /3 a; this is most adequate 
for MUB if we are mindful of Laplace’s advice to assign equal 
probabilities to alternatives about which we have no informa¬ 
tion EO). Second, after measuring all d -I- 1 bases, we know 
the density operator completely. 


d+l 

P = ^ Pa - 1 = 



alternatively, we write 


(2.5) 


^ d+l d-l 

^ctk hlafc ; with HJak Pak 

a—1k—0 


1 

d' 

(2.6) 
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While the mappings g i —pak and pak > Qa are linear, the 
mappings p^k ^ Q va Eqs. •HI or ( |2.6| l appear to be affine. 
In fact, these are linear mappings too, since we can make use 
of identities such as 


d+l d-1 




d-1 


a.k 1 


OL—l j,k—Q 


'^a.k — ^ ^ Pc^j I ^jk ^ j 5 


(2.7) 


if we wish. 

Here is a caveat. The linear map pak Qa in ( |2.51 l for a 
single basis continues to yield physical estimators if one re¬ 
places the probabilities pak by the corresponding relative fre¬ 
quencies from an experiment with a finite sample. The result¬ 
ing single-basis estimators have no statistical bias. For the full 
tomography, however, one cannot replace the probabilities in 
( |2.6| l by relative frequencies. If one does, one gets mock esti¬ 
mators that have no statistical bias but are unphysical, whereas 
physical estimators are statistically biased since they require 
a suitable nonlinear mapping from the relative frequencies to 
estimated probabilities. See Refs. lEiiia for various aspects 
of these matters. 

The second version of Eq. ( |2.5| ) exploits a well-known im¬ 
portant geometrical property of density operators: the differ¬ 
ence between any g and the completely mixed state ^ 1 is a 
traceless Hermitian operator, and these operators are elements 
of a — 1 real vector space. The sum over a in Eq. or 
( |2.6| > is a sum over hyperplanes in this vector space; there are 
d + 1 of them, each {d — l)-dimensional since w^k = 0- 
The hyperplanes are orthogonal in the sense of 


d-1 


i-1 


Tr Wak^ak X! ^/3/n/3Z = 0 , foT a ^ P , (2.8) 


\k^0 


1^0 


and in 


d+l d-1 


Tr(pp') = ;i + y]y] 


WakW^k : 


ct—l k—0 
d+ld-1 


Tr[(p-p')'] 

a—1 k—0 

we recognize the Euclidean metric of the vector space. 


(2.9) 


C. Linear inversion for incomplete MUB tomography 

The single-basis estimator in Eq. (in is obtained from 
the right-hand side of Eq. (in by setting P /3 i—^ for 
the basis indices /3 that are different from the privileged in¬ 
dex a. Now, if an intermediate number of bases from the 
set of MUB has been measured—say, the first M bases with 
1 < M < d + 1—the unbiasedness of the bases and the lin¬ 
earity just noted suggest to set i—> ^ 1 for the d + 1 — M 


unmeasured bases, which yields the Mth estimator of unbi¬ 
ased linear inversion (ULIN), 


a—l ^ ^ a—1 

^ M d-1 

— ~d^ ^ ^ ^ ^ ^ ^ak • 
a—1k—0 


M - 1 

( 2 . 10 ) 


We can also regard the truncation of the a summation as ef¬ 
fected by a replacement of the coefficients Wak in Eq. ( |2.6| ) in 
accordance with 


Wak 


Pak - ^ for 1 < Of < M , 
0 for a > M . 


( 2 . 11 ) 


Yet, irrespective of how we regard the map of the mea¬ 
sured probabilities pak with a < M onto the ULIN estimator 
the map simply amounts to estimating the probabilities 
of the unmeasured bases in the unbiased manner—faithful to 
Laplace’s advice, so to say: 

ppi = Tr(n^, ^ for /? > M . (2.12) 

All other estimators are characterized by the non-zero opera¬ 
tor they assign to the unmeasured part in 

1 v^) 

e-^l) , (2.13) 

/ unmeasured 

where 

^ \ (M) d-l-1 d-1 

p—= E E (2-14) 

/ unmeasured ^—M-\-l l—O 


^M) ^M) , 

^'other “ + 


is the estimator for the contribution of the unmeasured bases 
to the sums in Eqs. ( |0] l and ( |Z6l ). 

Before proceeding, six comments are in order. First, the 
map g i—>■ projects g — onto the span of the first M 

hyperplanes in the vector space discussed above; applying the 
map a second time has no effect. This projection property 
implies 


VeULlN t' / 


1 

d 


M d-1 

+ E E : 

ct—1 k—0 


(2.15) 


for any two density operators g and g'. 

Second, we note that, since Eq. ( |2.11[ ) involves only the 
measured bases, it is also well defined for systems where M 
MUB exist, but this set cannot be extended to d -f 1 bases. 
This is the case, e.g., for d = 6 and M = 3 where only sets of 
at most three MUB have been found thus far Il43lj48l . There 
are also non-extendable sets for prime-power dimensions Eol, 
although ULIN estimators for these sets suffer from the same 
problems as the ULIN estimators for subsets of full sets (see 
Sec. m- they are often unphysical. 
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Also, whether the unmeasured bases are pairwise unbiased 
among themselves, or unbiased with some of the measured 
bases, is not important for the maximum property of the fol¬ 
lowing third comment, nor for the physical least-bias estima¬ 
tors of Secs. |IVB] and [V^ If one wishes, one can choose the 
unmeasured bases such that they have largest average distance 
(see Ref. P9l ), so that they resemble a set of MUB as best as 
they can. 

Third, it is also worth noting that, of all estimators consis¬ 
tent with the measured probabilities, the ULIN estimator 
maximizes the Shannon entropy 

d+l d-1 

H^^\q) = - ^ (2.16) 

P=M+l 1=0 

associated with the unmeasured probabilities, which must be 
calculated according to Eq. ( |2.3| l. The ULIN estimator has this 
property by construction, since it gives the uniform estimates 
of Eq. p.l2| l for the ppis. 

For simplicity, we are taking some liberties with the Shan¬ 
non entropy in ( |2.16| i: the probabilities are normalized to unit 
sum for each /3, rather than to total unit sum and we are using 
the natural rather than the binary logarithm. This betrayal of 
the pure doctrine is of no consequence, however. 

Fourth, irrespective of the unbiased way of estimating the 
ppis, the ULIN estimator has a statistical bias; recall the per¬ 
tinent remarks in the Introduction. 

Fifth, as we shall see, after measuring M bases, we ac¬ 
tually have some information about the d + 1 — M unmea¬ 
sured bases. We are not really faithful to Laplace’s advice as 
long as we do not account for this information properly. For 
M = 1 we recover the estimator ( |2.4| i with a = 1, whereas 
= p is the full-tomography case of M = d + 1 in 
Eq. The reliability of for these limiting values 

of M is, however, misleading, and so is the self-suggesting 
probability estimation p.l2| l. As we shall see—and this is 
one main objective of this work—the ULIN estimator is not 
assuredly physical for intermediate M values, except in the 
qubit situation of d = 2. That is, ^ 0 is possible. The 
correct implementation of Laplace’s advice is not achieved by 
the simple-minded estimates in ( | 2 . 12 | i. 

Sixth, there is no other linear map pak ^ that could 

be used instead of p^k '—>■ ^lin if require for consistency, 
as we must, that pak = Tr(^^^ for the M bases that 
have been measured. There are, of course, legitimate state es¬ 
timators that are consistent with the known probabilities 
for bases a = 1,2,..., M —namely the physical “other” es¬ 
timators of Eq. p.l3| l —and these estimators make up a con¬ 
vex set; correspondingly, there is a convex set of estimated 
probabilities ppi for (3 = M + 1, ... ,d + 1. The ULIN esti¬ 
mator of Eq. ( |2.10| l and the unbiasedly estimated probabilities 
Ppi = 1/d may or may not be in the respective sets. When 
they are outside, we have to choose the best state estimator in 
accordance with a suitable optimality criterion. More about 
this in Sec. HY] 


III. ASPECTS OF INCOMPLETE MUB TOMOGRAPHY 
A. ULIN estimators for a single qubit (d = 2) 

As a first illustration, let us take the elementary example of 
a single qubit (d = 2). The three MUB here are just the eigen¬ 
states of the standard Pauli operators <Jy, and cr^, so that 
IIio = ^{t+a^) = na,+,..., Hal = 1(1 - (T^) = and 

Q = ^LIN =2^^+ SxCTx + Syay + s^a ^), (3.1) 

where s = (sx, Sy,Sz) is the three-dimensional Bloch vector 
with, e.g., Sx = Trigax) = ^iPx+ - Px-)- The eigenvalues 
of g are |(1 ± |s|), and the von Neumann entropy 

S'(p) = — Tr(pln p) (3.2) 

grows monotonously with shrinking Bloch-vector length |s|. 

When a single basis is measured (M = 1), say that of ax, 
only Sx is known, and we have 

d^JUN = ^(1 ■ (3.3) 

Clearly, this physical estimator maximizes the von Neumann 
entropy, since py± = pz± = | give the shortest Bloch vec¬ 
tor consistent with the known probabilities px± = |(1 ± Sa,). 
Likewise, the ULIN estimator for two measured bases, 

/o') 1 

PuLiN 2 ^^ “ 1 “ ^y^y) ; (3-4) 

which gives the estimated probabilities p^^ = Pz_ = 5 for 
the unmeasured CTz, maximizes the von Neumann entropy 
regardless of the values for the observed probabilities px^ 
and pyj^. 

These properties of the ULIN estimators for a qubit are 
consequences of the simple spherical geometry of the Bloch 
ball and the orthogonality of the planes corresponding to the 
distinct single-qubit MU Pauli observables. Such lines of ar¬ 
gument were employed for justifying the optimality of MUB 
tomography in higher dimensions when picturing the 

results obtained with a finite number of copies using “fuzzy” 
hyperplanes, their mutual orthogonality makes the uncertainty 
hypervolume particularly small. 

So, for a qubit, the unbiased unmeasured probabilities of 
the ULIN estimator, as quantified by the Shannon entropy, op¬ 
timize (maximize) the von Neumann entropy, much like the 
unbiased MU observables optimize (minimize) statistical un¬ 
certainty. These two properties complement each other. In 
fact, for any d > 2 and M, if g is equal to a projector IVak 
of one of the measured bases, then the ULIN estimator is the 
state itself according to the defining relation In view of 
all this evidence, one might expect this mutual compatibility 
to extend to higher dimensions—and so arrive at the ULIN 
estimator of Eq. ( |2.10[ ). However, we shall see that this some¬ 
what naive approach fails for d > 2 . 
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B. ULIN estimators for a single qutrit (d = 3 ) 


For the d = 3 case of a qutrit, we use this set of four MUB; 




(3.5) 




^1 0 0 ) 
04 = I 0 1 0 
.0 0 1 , 


where the columns are the probability amplitudes of the basis 
kets with reference to the fourth basis, and q = is the 

basic cubic root of unity. The 3x3 matrices that represent the 
respective single-basis estimators of Eq. (|2.4|i are 


1 


1 zi 

pi = Wzi 1 ^1 

^ \zi z; 1 


1 


1 


P 2 = g 1 ^5 


Z 2 q^zl^ 
1 q^Z2 


.qz2 qz2 


P3 = 


1 


23 92n 

J 

(P40 

0 


1 qz3 

: P4 = 

0 

P41 

0 

q^zi 1 

\ 

^0 

0 

P 42 J 


(3.6) 


where 


Za='^ q^Pak = Tr(p Za) , with Za^'^ q'^Bak ■ 

k—0 k—0 

(3.7) 

We have Pak = ^ + Wak = 5(1 + q~''za + q'^z’^) for the 
probabilities in terms of the z^s, and the unitary 3x3 ma¬ 
trices for the pairwise complementary observables Z^ are 

/O 0 

= 1 0 0 I , a = 1,2,3, 

\0 0 


'10 0 
Z 4 = \ 0 q 0 
,0 0 q^ 


(3.8) 


We note that pa = |(ll -f z* Zq, -f ZqZ^), and the replace¬ 
ment H-7 1 for an unmeasured basis is here simply im¬ 
plemented by zp I—)■ 0. 

Consider now a qutrit state for which the density operator 
is an incoherent mixture of two projectors, one each from the 
first and the second bases. 


g — Aillij + A2n2fe , 


(3.9) 


with 0 < Ai = 1 — A 2 < 1. As an immediate consequence of 
Trjllafc Zp) = Sapq^, we have here zi = Aiq^, Z 2 = X 2 q'^, 

and Z 3 = Z 4 = 0. It follows that = ^uun = ?{jun = Q, 
so that all ULIN estimators are genuine density operators for 
this particular g. 

Matters are quite different for 


for which zi = Z 2 = 2:3 = and Z 4 = — 59 ^; all gaS have 

eigenvalues 0, While = 0i and = Q are den¬ 

sity operators, as they are always, here ^un and ^lin are 
not acceptable estimators; they have negative determinants of 
— ^ and -“Y§g, respectively. 


The last example illustrates the statement in Sec. IIC that 
“the ULIN estimator is not assuredly physical for intermediate 
M values,” here for 1 < M < d + 1 = i. Analogous exam¬ 
ples can be constructed for higher dimensions; see Sec. [Ill C 


Regarding the other assertion in Sec. |IIC| — that there is no 
consistent “linear map p^k ^in ^ that could be used instead 
of Pak '—>> —we shall now argue that one gets a contra¬ 

diction when assuming that there is such a map. Lor this pur¬ 
pose, we look specihcally at the case of M = d = 3, but the 
argument can clearly be modihed for other dimensions d > 2 
and other intermediate M values. 

Suppose we have measured the probabilities of bases 1, 2, 
and 3 but lack data for basis 4. Then we know the values 
of zi, Z 2 , and Z 3 , which determine all off-diagonal elements 
in and we need a linear (or affine) map (zi, Z 2 , Z 3 ) 1 —>■ 


d3) 


for the estimation of the diagonal matrix elements. As an 


immediate consequence of the linearity, we note that convex 
sums of (zi, 22 , 2 : 3 ) values yield respective convex sums of 


43) 


values, 

zf\\zi + A'zj, Az 2 -I- A'z 2 , Az 3 -I- A'zg) 
= Xz^^\zi, Z 2 , Z 3 ) + A'zf^(zj, 22 , 23 ), 


(3.11) 


for 0<A=1 — A'<1. Lurther, for 21 = 22 = 23 = u with 
0 < |u| < we have 


u 0> 

= I w* • 0 
0 0 

= A I • 0 

0 0 


(3.12) 


+ A' 


0 ^ 


with = u/\u\ and A = | -f |u|, and the dots stand in for 
the yet-unknown diagonal entries. In the convex sum, both 
matrices must have p 4 o = P 4 i = \ and P 42 = 0 on the diago- 

nal, which implies the same values for this ^ip). We conclude 
that 


zf\u,u,u) = ^(1 + 9 ) = ■ 


(3.13) 


which, by convexity, also holds for it = 0 = ^it' + 3(—it'). 
Analogous arguments for 


zi = 922 = q^Z4 = It : ^^^=(000 

.It 0 3 


0 It* 


(3.14) 


and 


1 

-1 

0 \ 



-( 3 ) 
Plin — 

r 

0 

o\ 

-1 

1 

0 , 

( 3 . 10 ) 

Zi = 9^22 =qz 3 = u : 

0 

1 

2 

11 

VO 

0 

0 / 




VO 

u 

y 


( 3 . 15 ) 
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establish 


2,3,...,d+l, the nonzero matrix elements of 


z^i\u,q^u,qu) = + q^) = - ]^q , 

zf\u,qu,q^u) = + = _i . (3.15) 

Obviously, there is a contradiction for zi = Z 2 = = 0, as 

we cannot have —^ and —\q and also —\q^ for 24 ^^( 0 , 0 , 0 ). 
Case closed. 


Qa 

a 


'WaO 0 
0 Wal 


0 

0 


6^' 


L 0 


d—1- 


= Bi 


"0 * * • • •■ 

* 0 * • • • 

!t! * 0 • • • 




(3.20) 


C. ULIN estimators for higher dimensions 


The d = 3 example of Eq. ( |3.10| l is easily generalized to 
any dimension d. The ULIN estimator for M = d, 


-id) _ 

PuLIN — 


/ 

'- 1/2 

0 

0 


V 


- 1/2 0 0 ...\ 

1/d 0 0 ... 

0 1 /d 0 ... 

0 0 1 /d 


^0, 


■■■/ 


(3.17) 


is unphysical for all d > 2 because it has the negative eigen¬ 
value 1 /d — 1 / 2 . 

Yet another example of unphysical ULIN estimators for 
l<M<d+lis provided by a projector to a superposi¬ 
tion of two states from the hrst basis [the gs of Eqs. ( |3.10| l 
and p.l7| l project on superpositions of two states from the 
(d + l)th basis]; that is, 


Q = (l^/io) + IV'ii))^((^io| + (V^iil) 


= Bi 


- 1/2 1/2 0 •••■ 
1/2 1/2 0 ••• 
0 0 0 ••• 


B 


t 

1 f 


(3.18) 


where Si = ( {ipio) IV'ii) IV'id-i) ) is the row of kets 
from the hrst basis, as in Eq. (|33), and S| is the adjoint col¬ 
umn of bras. Note that we are using square parentheses to 
denote matrices expressed with the hrst basis, as opposed to 
the usual computational basis. 

Here, 


01 — ^LIN — Bl 


1/2 

0 

0 


0 0 
1/2 0 
0 0 


bI (3.19) 


accounts for all the diagonal elements whereas, for a = 


Here, the matrix to the ath basis is diagonal with the fcth en¬ 
try equal to Wak = Re((V’afc|V'io)/(V'afe|'0ii)). whereas the 
matrix to the hrst basis has null entries on the diagonal and 
some, if not all, off-diagonal entries are nonzero (indicated by 
the symbol *). Then, for M = 2,3,..., d, the d x d matrix in 


i?ULIN — dJl 


1/2 


* * • • 

1/2 * 

* 0 • • 


B\ 


(3.21) 


* * * • • • 0 . 


has Hermitian 2x2 submatrices 


* * 
* 0 


^ 0 with nonzero 


off-diagonal elements and one or two vanishing diagonal el¬ 
ements. Such a 2 X 2 matrix has a negative determinant and 


cannot be positive semidehnite. It follows that ^ 0 is 
unphysical for intermediate M values. 

Lor illustration, we take once more the qutrit MUB in 
Eq. (|^, for which 


^ 1 ) _ 

PuLIN — 


1 

2 


1 0 0 
0 1 0 
0 0 0 


^2) _ i 

PuLIN — « 
6 


3 1 q 
1 3 g , 
q^ q^ 0 


-(3) 1 

PuLIN — n 

6 

■3 2 - 1 ' 

2 3-1 

^4) 1 

7 PULIN — 0 

'l 1 o' 
1 1 0 

-1 -1 0 

Z 

0 0 0 


(3.22) 


Indeed, ^lin ^nd are unphysical. 

In summary, then, we have seen in so many examples that 
ULIN estimators can be unphysical. The unbiased estima¬ 
tors Qp = 2 ^ are not permissible for /3 > M and a nonzero 
“unmeasured” term is needed in Eq. ( |2.13| ) for a physical esti¬ 
mator In this situation, the Md measured probabilities 

Pak with a < M impose constraints on the (d -|- 1 — M)d un¬ 
measured probabilities, which exclude the use of ppi = 7 for 
(3 > M. It follows that we have some, if limited, knowledge 
about the observables associated with the unmeasured bases: 
our best guess is not that all outcomes are equally likely for the 
future von Neumann test of a basis not measured as yet. Put 
differently, physical state estimators for 1 < M < d -f 1 typi¬ 
cally possess some bias in the unmeasured probabilities, such 
that they deviate from the uninformative uniform distribution. 
This is the reason why we cannot implement Laplace’s advice 
by the naive estimates of Eq. (|2.12]i: it is simply not true that 




































7 


we have no information at all about the unmeasured probabil¬ 
ities. 

Harking back to the qubit estimators in Eqs. CT and ( [I4l i, 
we observe that these also contain information about the un¬ 
measured bases because < 1 must hold for all 

physical estimators, so a known value for x restricts both y 
and z, and known values for x and y bound the acceptable 
z values. Since y = z = 0orz = 0, respectively, are always 
permissible estimates, the intermediate ULIN estimators work 
for the qubit. This feature, however, has no analog for d > 2. 

D. Nonpositivity of ULIN estimators 

Before moving on to discussing the proper choice for the 
“unmeasured” contribution to in Eq. ( |2.13| ), let us close 
the subject of negative eigenvalues of intermediate ULIN es¬ 
timators by a more quantitative study. We ask this question: 
what is the most negative eigenvalue Amin that pos¬ 

sibly have? 

Eor given p, the smallest eigenvalue of can be found 
by minimizing Tr(^f^N o') over all density operators a, and 
further minimization over p establishes 

Amin = minTr(^uN o’) = minTr(po=^,uN) > , • 

g,(T Q,(7 a 

(3.23) 

Here, the roles of p and cr are interchangeable thanks to the 
symmetry noted in Eq. ( |2.15| l, and the lower bound follows 
from Eq. ( |2.10| i in conjunction with Tr(pQ, o') > 0 for all a. 

The equality holds in Eq. ( |3.23[ ) only if Pa is rank- 

deficient, and this is the case for certain (d, M) pairs but not 
for others; see Eig. [ll Eor example, Amin = — \ can be es¬ 
tablished for {d,Mj= (4,2) by a simple explicit example. 
Whether tight or not, the lower bound in Eq. ( |3.23| l tells us that 
partial MUB tomography with large d and small M is very 
nearly unbiased. Eor example, if only two of the MUB are 
measured, the lowest minimum eigenvalue of the ULIN esti¬ 
mator cannot be smaller than —1/d, and would approach zero 
for large dimensions. Yet, even a small negative eigenvalue 
renders unphysical if this eigenvalue is different from 
zero within the numerical accuracy with which it is known. 

The extremal density operators p and cr minimizing the 
traces in Eq. ( |3.23[ ) are such that a projects on the smallest 
eigenvalue of p^^ and p projects on the smallest eigenvalue 
of if 0^)2 ^ so that we have the pair of equations 

^unO- = Tr(^uu^cr)cr, = Tr(CT^))2p)p. (3.24) 

We solve them by iteration: 

SI: Eor the current p, diagonalize set cr to the pro¬ 

jector onto the smallest eigenvalue. 

S2: Diagonalize update p by setting it to the pro¬ 

jector onto the smallest eigenvalue. 

S3: Repeat steps S1 and S2 until ( |3.24| l holds with the desired 
numerical accuracy. 



FIG. 1. The most negative possible eigenvalue Amin of a ULIN 
estimator in partial MUB tomography, for prime-power di¬ 

mensions in the range 3 < d < 13. The successive Amin values of 
M — 1, 2,..., d + 1 are connected by straight-line segments that 
guide the eye. The curves are shifted by —0.2 relative to the preced¬ 
ing d value. The vertical-axis labels are shown for d = 3 and d = 9 
on the left, and for d = 5 and d = 13 on the right. Each plot point 
represents the optimal value of Eq. l |3.23^ over one hundred repeated 
numerical searches for solutions of the equation pair ( |3.24^ with the 
algorithm discussed in the text. Filled circles (•) indicate values with 
“=” in Eq. l |3.23| l, including M — 1 always; empty circles (o) mark 
values with “>”, including M = d -|- 1 always. Eor M = d, where 
“>” is the case, we have the negative eigenvalue ^ — 7 of the exam¬ 
ple in Eq. l |3.17^ , shown by filled squares (■). 

Since we are minimizing a convex function over a convex 
domain, the possibility of landing at a suboptimal point on the 
boundary cannot be ruled out. Restarting the numerical search 
several times with different initial states is necessary. 

Figure shows these most negative eigenvalues for 
the prime-power Hilbert-space dimensions from d = 3 to 
d = 13. For d = 7 and M = 2, for instance, we find that 
Amin ~ —0.1394 while the example of Eq. ( |3.18| ) gives a most 
negative eigenvalue of about —0.1250. Further, we observe 
the following features: 

FI: The most negative eigenvalue is obtained for M = d, and 
is equal to ^ the negative eigenvalue of ^un in 
Eq. ( |3.17| l. These values are marked by filled squares 
in Fig. m 

F2: For 2 < M < d, there are always more M values for 
which “>” applies in Eq. ( |3.23| l than those for which “=” 
is the case. These values are marked by empty or filled 
circles, respectively. 

F3: Some deviations from a monotonic decrease of Amin as 
M increases from 1 to d are observed for some dimen¬ 
sions, notwithstanding the obvious trend. 

We leave it as a moot point whether or not these features are 
also present in dimensions higher than d = 13. Should feature 
El be generally true, then 

(3.25) 

would sharpen the inequality in Eq. (|3.23|l. Currently, this is 
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just a conjecture suggested by the evidence presented in Fig.[T] 
Whether it holds for all prime-power dimensions d, is perhaps 
of some interest for those who study the properties of MUB. It 
is, however, of no consequence for quantum state estimation, 
where the fact that Amin is negative for all intermediate M 
values matters, while the precise value does not. 


IV. LEAST-BIAS MUB INFERENCE 

A. Physical unbiased estimators and von Neumann entropy 

As we know, if the ULIN estimator is a bona fide 
density operator, it maximizes the Shannon entropy of the un¬ 
measured probabilities in Eq. ( |2.16| l. For d > 2, only the es¬ 
timators for M = 1 and M = d + 1 will surely also yield 
the largest von Neumann entropy of Eq. ( |3.21 l. The permis¬ 
sible ULIN estimators for \ <M<d+ldQ not generally 
maximize the von Neumann entropy over the convex set of 
the physical of Eq. (|2.131i: they are not the estimators 
of Ref. mi. 

A simple counterexample for d = 3 and M = 2 is enough 
to support this statement. We return to Eq. p.9| l and con¬ 
sider the equal-weight case of Ai = A 2 = so that 21 = ig-i, 
^2 = 5 ?^, and Z 3 = Z 4 = 0. The rank-2 operators g = 
^UUN = diuN = have eigenvalues 0, |(3 ± ^3) and 

von Neumann entropy 5 '(^lin) = 0.5157. The largest value 
of the von Neumann entropy is 0.6370, however; we obtain it 
for the with = -0.09466. 

It is, of course, hardly surprising that the Shannon and 
the von Neumann entropy are usually maximized by differ- 
ent ^2ier^- These are two different figures of merit, which 
serve different purposes; if we intend to measure the remain¬ 
ing bases of the set of MUB, the Shannon entropy is a useful 
quantity; by contrast, if we want to have a more universally 
applicable estimator, then maximizing the von Neumann en¬ 
tropy is a time-honored approach in Jaynes’s spirit lEHSa. 
That the two procedures coincide in the qubit case is nothing 
more than a result of the geometrical simplicity of the prob¬ 
lem. 


B. Physical estimators with least bias 

Yet, in the context of incomplete MUB tomography, it is 
clearly more natural to choose the “unmeasured” contribution 
in Eq. ( |2.13| l such that the Shannon entropy is maximized 
rather than the von Neumann entropy. For, the maximiza¬ 
tion of the Shannon entropy ensures that we accept 
the estimator whenever that is permissible—that is, whenever 
^LJN ^ 0—when that is not the case, we stay as close to 
naive unbiasedness as possible (more about this closeness in 
Sec.[^. In other words, it is not necessary to completely dis¬ 
card the notion of unbiasedness; a certain refinement is called 
for: we opt for the least-bias estimator which is the 


physical with the largest Shannon entropy for the un¬ 
measured probabilities, 

max ■ (4.1) 

e‘r>o 

We repeat, perhaps unnecessarily, that the least-bias estima¬ 
tor is different from the ULIN estimator only if the latter is 
unphysical. 

C. An efficient algorithm 

The numerical search for requires the maximization 
of H^^\q) over the convex set of physical 
over the set of all density operators g with = Tr(pnQfc) 
for a = 1,2,..., M. Other than that, g is only constrained 
by the defining properties of a density operator: p > 0 and 
Tr p = 1. This optimization problem is very similar to the one 
solved in Ref. m, where one maximizes the von Neumann 
entropy over the convex set of maximum-likelihood estima¬ 
tors. Indeed, the algorithm of Ref. nsi can be modified so 
that it applies to the current problem of finding ■ 

The figure of merit is the following function of p: 

M d-1 

Df,{g) = X!7’afcln[Tr(pnafc)] 
a—1k—0 

rf+1 d-1 

E E Tr(pnp/) ln[Tr(pn^;)] . (4.2) 

/3=M+1 1=0 

The first term is a mock log-likelihood for relative frequen¬ 
cies equal to the known probabilities p^k for a < M. The 
second term is the Shannon entropy I-l^\g) of multiplied by 
the nonnegative parameter p. 

Now, all physical maximize the p = 0 function 

Do{g), with 

M d-1 

maxUoCe) = = '^'^Pak Inpafe . (4.3) 

a—1fc—0 

Because of this extremal property of first- 

order contribution to D^{g) originates solely in the Shannon- 
entropy term, 

max U^(p) = ^) H-, (4.4) 

e 

where the ellipsis stands for contributions of order /i^ and 
higher. Therefore, the numerical maximization of D^{g) for 
sufficiently small p yields the least-bias estimator . Prac¬ 
tical experience shows that p « 10“^ is large enough for a 
noticeable difference between D^{g) and Do{g) and small 
enough to ensure that the density operator that maximizes 
D^{g) approximates very well. Usually, it is a good idea 
to numerically select a p that minimizes the fluctuations in 
the extremal p—from the iteration described below—to some 
pre-chosen precision. 
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The first-order response of D^{q) to a variation q i—> p + 6p 
is 


bD^[Q)=TT[bQW{Q)], (4.5) 


with 


M d-l 


w { q ) = Y.11 


Pock 


n 


Oik 


a—1k—0 

d+1 d-l 


Tr(£»na/e) 

d-l 

E En.. ln[Tr(pn^i)]. (4.6) 

/3=M+1 1=0 

Hence, the positivity constraint for g implies the following 


extremal equation for ; 

. (4.7) 

Since D^{q) has no local maxima, we can find by 

steepest ascent (“follow the gradient uphill”). The W (g) iden¬ 
tifies the gradient in the sense that 

bg(xW(g)g+gW(g)-2gTrlgW(g)] (4.8) 

gives the largest first-order change in D^{g). We ensure a 
positive increment of D^{g) by the iteration in accordance 
with 


Qn i'ra-l-l — 


{i+e{W{gn) - Tr[iy(e„)g„]ll}) g„(l + e{H^(g„) - Tr[iy(g„)e„]]l}) 
1 + e2 Tr({VF(p„) - Tr[iy(p„)p„]l}"p„) 


(4.9) 


where e is a small positive step size. To obtain the extremal 
solution one can start with the maximally-mixed state 

gn^i = 2 ^ and iterate Eq. ( |4.9| l for a small fi until the ex¬ 
tremal equations ( |4.7| i are obeyed to a satisfactory accuracy. 

We note, but do not elaborate this point, that one can speed 
up the convergence substantially by employing conjugate gra¬ 
dients ll54ll : each iteration step is then more costly (in CPU 
time) but that is more than compensated for by the much 
smaller number of steps. Further we note that, in cases where 
the iteration algorithm proceeds too slowly due to the com¬ 
plexity of the optimization problem for large Hilbert-space di¬ 
mensions, the sum of the two entropic functions in Eq. ( |4.2| l 
can alternatively be optimized with gradient-free methods, 
such as the Nelder-Mead “amoeba” method ll55]l or simulated 
annealing ll5^ . 


D. Examples 

For explicit examples, we look at a class of qutrit states 
comprising statistical mixtures of g as represented in p.lO[ ) 
with the maximally-mixed state for various admixtures 0 < 
w < 1: 

g>-^ gw = {'^-'w)g + w-. (4.10) 

As always, for M = 1 and M = 4, the ULIN estimator is 
always positive semidefinite. Beyond w = 0.2679 and w = ^ 
respectively for M = 2 and M = 3, the ULIN estimator is 
precisely the physical least-bias estimator that maximizes the 
Shannon entropy since the true state would be highly mixed. 
For other w values, the ULIN estimator possesses at least one 
negative eigenvalue. The relevant parameters for the unphysi- 
cal ULIN and physical least-bias estimators for two exempli¬ 
fying values of w are listed in Table |I] 


TABLE I. Comparison between the ULIN and least-bias estimators 
that maximize the Shannon entropy for Qw When p™ is the pure 
state defined by Eq. l |3.10t (w = 0), the least-bias estimator is the 
true state g. Also, note that for zi = Z 2 = Z 3 = —(1 — w)/2, 
Z4 = —(1 — w)q^/2. 



w 

= 0.1 

w = 0.2 


Zi 

23 

24 

Qvj 

-0.450 

0.225 -hi0.389 

-0.400 

0.200-hi 0.346 

^2) 

^ULIN 

0 

0 

0 

0 


-0.313 

0.156-hi 0.271 

-0.126 

0.063-hi 0.109 

^3) 

^ULIN 

-0.450 

0 

-0.400 

0 


-0.450 

0.174-hi 0.303 

-0.400 

0.100-hi 0.173 


V. COMPARISON WITH OTHER ESTIMATORS 
A. Estimators of the least-bias kind 

The Shannon entropy quantifies our ignorance about the 
outcomes of future projective measurements in the yet- 
unmeasured bases with P = M + 1, M + 2,... ,d + 1 from 
the set of MUB. Rather than ignorance we can equiva¬ 
lently quantify knowledge by the predictability of the fu¬ 
ture measurements. For an experiment with probabilities 
Po,Pi, ■ ■ ■ ,Pd-i 7 the entropic measure of predictability is 

d-l 

-Pent(p ) logd(pid), 


(5.1) 
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23 times magnification 



FIG. 2. A prototypical showcase of three types of physical least- 
bias estimators relative to the ULIN estimator for different qutrit true 
states (M = 3). The black triangle is the regular simplex for «4 in 
the complex plane, with its corners at Z 4 = 1, 5 , and (f. The blue 
curve borders the regions of permissible values for 04 . The examples 
featured here are for w = \ and state 1 of Table[n| The magnified 
view is for state 1 , where the three least-bias estimators are different; 
by contrast, all three are the same for ■ 


where the symbol p. stands for all the d probabilities. We have 
the extreme values 


Pent{p) 


0 

1 


ifPo =^1 = ••• =Pd-i = 
if Pi = Sii for a certain I, 


(5.2) 


and 0 < PentiP-) < 1 for all other sets of probabilities, which 
are among the defining properties of all permissible measures 
of predictability ll57l l58l . Convexity 

P{Xp. + X'p') < XP{p.) + X'P{p'.) (5.3) 


for 0 < A = 1 — A' < 1 is another important property of all 
predictability measures. 

In terms of Pent(p.), the Shannon entropy is 




d+1 


d+l-M- Y. Pent(Tr(pn^.)) 
P=M+l 


Ind, 


I that we could replace the figure of merit in Eq. (|4.2[) by 


(5.4) 



FIG. 3. A comparison of the lines of constant “distance” from the 
ULIN estimator for the three different predictability functions that 
are used as figure of merit for the least-bias estimators in Figs. 2 and 
4. We have P = 0 at the center of the triangle and P = 1 in each 
comer, and the contour lines indicate predictability values P = 0.2, 
0.4, 0.6, and 0.8, respectively. The green circles are for the purity 
measure of Eq. ( |5.6| l; the orange distorted circles are for the Shannon- 
entropy measure of Eq . l|5.1[ > and the purple hexagons are for the 
betting measure of Eq. l|5.7|l. 


Upon replacing Pent( ) in ( |53| ) by Ppur(), the that maxi¬ 
mizes D^{g) is simply another least-bias estimator than 
the one for Pent( )■ The two least-bias estimators are differ¬ 
ent since they refer to different quantifications of the bias re¬ 
garding the probabilities for the not-yet-measured bases. But 


neither is better than the other; they serve different pur¬ 
poses. Any other predictability P{p.) yields a corresponding 
least-bias estimator. Our preference, in Sec. IV for the g^s^^ 
associated with Pent(p ) is mostly for historical reasons— 
history of the subject, that is. 

Yet another example is provided by the predictability func¬ 
tions that equal the expected gain for various strategies of bet¬ 
ting on the future outcome. There is, in particular, the “linear 
bet” of Sec. 2.3.1 in Ref. ll5^ . which amounts to 


Pbet{p ) = max{po,Pi,P2} - mm{po,pi,P2} (5.7) 


M d-l 

Df,{g) = ^ ^Pafcln[Tr(pn„fc)] 

a—1 k—0 
d+1 

-P Y ^"ent(Tr(pn;3.)), (5.5) 

P=M+1 

with a corresponding minor change in Eq. ( |4.6| l, namely 
/i I—)■ p/lnd. 

Other measures of predictability could be employed as well. 
Eor example, the purity-based predictability 


in the qutrit case. 

We can regard Ppui (p.) in ( |5.6| ) as the normalized squared 
Hilbert-Schmidt distance of po,pi,... ,pd-i from the uni¬ 
form distribution p; = i. Likewise, Pentip ) is the normalized 
Kullback-Leibler divergence of the actual probabilities from 
the uniform distribution. Analogous remarks apply to any 
other permissible predictability P{p.). It follows that, irre¬ 
spective of which P{p ) we choose, the resulting least-bias es¬ 
timator equals the ULIN estimator whenever ^ 0- And 
when ^ then—as a consequence of ( |5.3[ ) —each 

sits on the border of the convex set of physical what¬ 

ever border point is closest to in the sense specified by 
the choice of P{p.). 























11 


These matters are illustrated in Fig. which graphically 
shows the situation for qutrit MUB where M = 3 bases are 
measured. As it turns out, numerical experience shows that 
different types of least-bias estimators give estimators that are 
quite close to each other, hinting that the choice of measure 
for quantifying bias typically influences the resulting estima¬ 
tor only mildly as long as these nonpathological choices lead 
to the same extremal solution in the absence of the quantum 
positivity constraint. Such differences would not matter much 
in the presence of statistical fluctuations and other experimen¬ 
tal error sources. For various qutrit states, the two kinds of 
least-bias estimators introduced in this section are essentially 
identical. For completeness, in Fig.[^we present a comparison 
of the different figures of merit for the least-bias kind. 


B. Other estimators 

The least-bias estimators thus constructed are rank- 
deficient whenever ^ 0. It may, however, be desirable 
to use estimators that usually have full rank—they are more 
robust in the sense that a slight perturbation, or inaccuracy in 
determining them, does not render them unphysical. Other 
reasons for avoiding rank-deficient estimators have been put 
forward as well Il59l . 

From the plethora of physical we consider three 

choices: the estimator with the largest von Neumann en¬ 
tropy m, the Bayesian mean estimator for a uniform prior 
on the unmeasured probabilities 1 ^ . and the estimator hav¬ 
ing the largest minimal eigenvalue MM- 

These estimators simply differ in the unmeasured probabil¬ 
ities, as is the case for any estimator consistent with the prob¬ 
abilities of the already-measured bases. To streamline pre¬ 
sentations, Fig. and Table [H] illustrate the estimators for the 
M = 3 case of a qutrit. 


VI. CONCLUSIONS 

In summary, we have examined whether the ideal Laplace- 
type linear estimator that assigns equal and unbiased proba¬ 
bilities to all the outcomes of unmeasured bases from the set 
of MUB is a physical thing to do in incomplete MUB tomog¬ 
raphy. The answer is definite: such unbiased estimation does 
not work as a general recipe. 

As a natural adjustment to the original inference method 
that follows a blind application of Laplace’s notion of indif¬ 
ference, we recommend the use of the least-bias estimator. 
It gives the closest-to-uniform unmeasured probability distri¬ 
butions by maximizing the Shannon entropy with due atten¬ 
tion to the constraints imposed by the measured probabilities. 
We also supply a simple iterative algorithm for computing the 
least-bias estimator. 

We compared this least-bias estimator with alternative esti¬ 
mators that quantify the bias differently, and concluded that 


these alternatives are equally useful for practical purposes. 
Since all these least-bias estimators are rank-deficient when¬ 
ever thej^jiie different from the linear estimator of ideal 

☆ ULIN 


FIG. 4. Pictorial comparison of the ULIN estimator with three phys¬ 
ical full-rank estimators; the estimators with largest von Neumann 
entropy and largest minimal eigenvalue, and the Bayesian mean es¬ 
timator. The rank-deficient three least-bias estimators of Fig. 2 are 
also indicated; they sit on the blue border of the region of permissi¬ 
ble states and are too close to each other for telling them apart. As in 
Fig. 0 the examples are for qutrit states with M = 3, and the triangle 
is the simplex for 04 . The examples featured here are state 2 and 
form = j of Table [n| 

Laplace-type, we also took a look at three different estimators 
of full rank. As illustrated by examples, the full-rank estima¬ 
tors are indeed different from the least-bias estimators, but are 
certainly acceptable as consistent bona fide estimators. 

Armed with these insights, one could now study questions 
such as after measuring in M bases, which basis is optimal 
for the next von Neumann test? This matter and others are, 
however, beyond the scope of the present article. 

ACKNOWLEDGMENTS 

Many of the ideas in this paper originated at the Work¬ 
shop on Mathematical Methods of Quantum Tomography at 
Fields Institute (Toronto) in February 2013. Z. H., J. R., and 
Y. S. T. are grateful for the support of the European Social 
Fund and the state budget of the Czech Republic [project No. 
CZ. 1.07/2.3.00/30.0004 (POST-UP)], the Grant Agency of the 
Czech Republic (Grant No. 15-031945), the IGA Project of 
the Palacky University (Grant No. IGA PrF 2015-002), and 
the sustainability of Post-Doc positions at Palacky Univer¬ 
sity. L. L. S. S. acknowledges the support from UCM-Banco 
Santander Program (Grant GR3/14) and helpful discussions 
with Markus Grassl. H. K. N.’s, J. H. C.’s, and B.-G. E.’s 
work is funded by the Singapore Ministry of Education (partly 
through the Academic Research Eund Tier 3 MOE2012-T3- 
1-009) and the National Research Eoundation of Singapore. 
H. K. N is also funded by a Yale-NUS College start-up grant. 



O maximal Shannon entropy 
V minimal purity 
X minimal linear bet 

# maximal von Neumann entr( 
A largest minimum eigenvalue 

♦ Bayesian mean 



State 2 





12 


TABLE II. Values of 24 for physical estimators to three un¬ 
physical ULIN estimators for qutrit states with three measured 
bases (d = M = 3). These data are used in Figs. 1^ and 1^ The 
values of 21 , 22 , and 23 that specify the respective ^^s are 
(i) 2 i = 22 = 23 = —3/8 for (ii) zi = 0.160 — i 0.321, 

22 = 0.571 — i 0.192, and 23 = 0.314 + 10.165 for state 1; 

and (hi) 21 =—0.345 + i 0.0574, 22 = 0.303 + 10.328, and 

23 = 0.00057 — i 0.294 for state 2. The rows report the correspond 


ing three 24 values for the least-bias estimator of S ecs. 
and 


IV B 


IV C 


for the alternative least-bias estimators of Sec. 


VA 


^PUR 


and for the full-rank estimators with largest von Neumann en¬ 


tropy and largest minimal eigenvalue, and and for the 


Bayesian mean estimator . For information, the last row contains 
the 24 values of the actual states that were used to obtain the values 


of 21 , 22 , 23 for the 



24 

state 1 

24 

state 2 

24 


0.067 + 10.106 

0.080 + 10.299 

0.073 

-i0.136 


0.067 + 10.106 

0.093 + 10.295 

0.073 

-i0.136 

+3) 

t^BET 

0.067 + 10.106 

0.128 + 10.289 

0.080 

- i0.132 

-<3) 

5vN 

0.120 + 10.208 

0.090 + 10.309 

0.104 

- i 0.204 

-<3) 

^MINEIG 

0.187 + 10.325 

0.003 + 10.438 

0.122 

- i 0.283 

^3) 

^BM 

0.176 + 10.305 

0.021 + 10.418 

0.122 

-i 0.279 

actual state 

+i 0.324 

-0.146 + 10.601 

0.256 

-i 0.337 
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